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Abstract 



2 ■ We extend the standard theory of cosmological perturbations to homogeneous but anisotropic 

c/3 , universes. We present an exhaustive computation for the case of a Bianchi / model, with a 

^ ' residual isotropy between two spatial dimensions, which is undergoing complete isotropization 

at the onset of inflation; we also show how the computation can be further extended to more 
general backgrounds. In presence of a single inflaton field, there are three physical perturbations 
' (precisely as in the isotropic case), which are obtained (i) by removing gauge and nondynamical 

■ degrees of freedom, and (ii) by finding the combinations of the remaining modes in terms of 

which the quadratic action of the perturbations is canonical. The three perturbations, which 
^ later in the isotropic regime become a scalar mode and two tensor polarizations (gravitational 

Q ■ wave), are coupled to each other already at the linearized level during the anisotropic phase. 

, This generates nonvanishing correlations between different modes of the CMB anisotropies, 

^ ' {O'lm^i'm') ^ ^u'^mm'j which Can be particularly relevant at large scales (and, potentially, 

be related to the large scale anomalies in the WMAP data). As an example, we compute the 
^ . spectrum of the perturbations in this Bianchi / geometry, assuming that the inflaton is in a slow 

^ I roll regime also in the anisotropic phase. For this simple set-up, fixing the initial conditions for 

the perturbations appears more difficult than in the standard case, and additional assumptions 
seem to be needed to provide predictions for the CMB anisotropies. 



1 Introduction 

The study of Cosmic Microwave Background (CMB) full sky maps from the WMAP experiment [1] has led 
to some intriguing anomalies which seem to suggest that the assumption of statistical isotropy is broken on 
the largest angular scales [2j . These anomalies include an alignment of the moments in the lowest multipoles 
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dubbed the 'axis of evil', an asymmetry in the power between the northern and southern echptic hemispheres 
and an apparently non-Gaussian excursion in the southern galactic hemisphere known as the 'cold-spot'. 
Another well-known puzzle in the observations has been the lack of power in the quadrupole which had also 
been noted in the COBE-DMR maps [3J and is still present in the latest WMAP data [1]. 

The significance of the anomalies has been debated extensively in the literature (see e.g. [5]) with some 
reported effects more significant than others. The difficulty in quantifying exactly the importance of any 
effect is due to how correctly the a posteriori probability of observing it is estimated. It is clear however 
that the anomalies implying an overall anisotropy in the data are more significant than the lack of power 
in the CMB quadrupole. A natural explanation for the observed anomalies may be some form of as yet 
undetermined systematic or foreground signal which is not being taken into account properly in the data 
reduction producing the final maps [B]. However, a conclusive explanation along these lines has not been 
put forward yet. It is therefore legitimate to ask whether the observed anomalies may be an indication of a 
departure from the standard cosmological model. In order to find convincing evidence on such deviations, 
one should study the CMB properties predicted by any specific model, with a particular emphasis on those 
which go beyond the ones characterizing the standard picture. For instance, if the signal is statistically 
anisotropic, crucial information will be encoded in the off-diagonal correlations of the spherical harmonic 
modes with (a^m^^'m') ^ ^ii'^mm' ^ which are predicted to vanish in the standard case. 

Since, aside the anomalies [2j, CMB data strongly support the general theory of inflation, we focus here 
on departures from the standard picture which can however be reconciled with the inflationary framework. 
For example, the low quadrupole power has motivated a number of studies to explain the observed cut-off 
within the context of inflation O [HI [9l [10] (ref. [llj further extended the model [7] , to account for a spatial 
asymmetry of the CMB perturbations at large scales). These models assume a short period of inflation, 
leading to a fine tuning problem in the standard inflationary picture since the slow-roll regime responsible 
for the observed ~ 60 or so e-folds is an attractor in the inflaton phase space and thus the probability of 
inflating for a large number of e-folds is much larger than that of inflating for only 60. Recently, other 
groups have argued that inflation should be short-lived [13] and in particular models of inflation on the 
string theory landscape seem to predict a highly suppressed probability of large A^ [14J . Coupled to the fact 
that our observed universe requires A^ ~ 60, these arguments imply that inflation must have lasted 'just 
long enough'. Theoretical arguments aside, a limited amount of inflation would certainly lead to a richer 
phenomenology than a prolonged one. The reason why inflation was originally postulated is that it leads to 
the isotropic and homogeneous Friedmann-Robertson- Walker (FRW) universe starting from rather generic 
initial conditions. If inflation lasted many e-folds, then any trace of the pre-existing universe is inflated 
away to scales much larger than the present day horizon. On the other hand, a limited amount of inflation 
could have left some trace in the data, particularly at large scales, that would appear as anomalies within 
the standard inflationary picture. 

With these considerations in mind, in this work we perform an initial step towards the study of the 
CMB anisotropics taking into account a pre-existing phase before inflation, in which the universe is not yet 
of the FRW type. Motivated by the suggested alignment of the large scale multipoles along a common axis, 
we focus on models which are anisotropic at the onset of inflation. For simplicity, we still assume that the 
universe is homogeneous. As we will see, even this very minor departure from standard cosmology requires a 
signiflcant extension of the well consolidated framework for the computation of CMB anisotropics in a FRW 
universe. Homogeneous but anisotropic cosmologies have long ago been classified into equivalence classes 
known as the Bianchi types [15]. With the possible exception of type /A, Bianchi universes with a positive 
cosmological constant evolve towards an asymptotic (isotropic) de Sitter stage [16]. The situation is more 
complicated when the cosmological constant is replaced by an infiaton field (which one has to do, in order 
to have a finite amount of infiation). Also in this case, however, one finds that inflation takes place for a 

''See however other inflation models that naturaUy suppress the overall number of e-folds, e.g. [12) . 
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wide variety of circumstances, leading to the isotropization of the space. □ 

It was shown in [TH] that a late time contribution of the Bianchi Vllh form to the temperature 
anisotropics can account for several of the WMAP anomalies; more in general, Bianchi models with gen- 
eralized backgrounds driven by anisotropic stresses have been studied extensively in the literature (see e.g. 
|19|, I20j). These models require a particular mechanism actively driving the anisotropy, which can for this 
reason survive for a much longer period than in the case considered here (and, possibly, result in stronger 
observational effects) . In this work we consider only the case consistent with the simplest models of inflation 
where the background is initially anisotropic and is isotropized at the onset of inflation, [fl We expect to 
recover a standard power spectrum at the scales that leave the horizon after the isotropization has been 
achieved, but nonconventional results at larger scales. 

Thus, the anisotropy is restricted to the initial conditions, in the form of a power spectrum of the pri- 
mordial perturbations which depends also on the directionality of the modes, and not only on the magnitude 
of their momentum k (as in the isotropic case). On the contrary, the transfer functions, relating the late 
time temperature anisotropics to the primordial perturbations, are the standard ones, since the propagation 
of the modes, soon after the onset of inflation, occurs in an isotropic background. This simplifies greatly 
the calculation of the correlations induced in the CMB. H We present this computation in Section O in the 
case in which there is a residual anisotropy between two of the three spatial directions. 

Namely, we compute the correlation between different CMB modes, under the assumption that (i) 
the primordial perturbations are Gaussian, and (ii) their power spectrum depends on the magnitude |k| 
and on the angle between the momentum and the anisotropic direction. We do not impose any further 
assumption in this computation, so that it holds for any Bianchi model with a residual 2d isotropy and an 
asymptotic standard inflationary behaviour. As we anticipated a few paragraphs above, the main result is 
a non vanishing correlation between off diagonal modes, see eq. ([7j). The specific relations found between 
the different non vanishing correlators can in principle allow for a detailed data analysis within any given 
model. We do not perform this analysis in the present work (since, as we mentioned, our main result is 
the extension of the computation of the primordial perturbations to such non FRW geometries); however 
we include the relations ([7]) , since they offer a model oriented way of extracting information from the CMB 
data. This information is lost in the standard diagonal Ci statistic analysis. The observed correlations 
can be fit with model-independent templates (see e.g. [271 |l]) but these are not calculated a priori from 
any particular theory and thus the information contained in the CMB will not constrain any fundamental 
parameter. 

In the remainder of the work, we present the detailed computation of the primordial perturbations within 
a Bianchi / model with a residual 2d isotropy, which is the simplest deviation from the FRW geometry. The 
line element is of the form 



and late time isotropy (a = b) is achieved due to an inflaton scalar field cj) . In Section [3] we study the 
background evolution for this model, under the assumption that the infiaton is initially in a slow roll 
regime. As it is well known, the anisotropy is rapidly damped away, on a timescale corresponding to the 

^ A list of works discussing this issue, and several aspects of inflation in Bianchi models can be found in the review [17j . 

''We also assume that the only source is the inflaton field. Ref. [21] studied the background evolutions in cases in which an 
anti-symmetric tensor is also dynamically relevant. 

'^Ref. [22] studied the perturbations of a scalar field assuming that a small anisotropy (corresponding to the ratio h/H in 
our equation (|14p) is present during infiation, and that it then decays away at the end of infiation. Such perturbations are then 
converted into metric perturbations through the mechanism of modulated perturbations [53]. This work presents analytical 
results for the perturbations using an expansion series in the anisotropy parameter, which, in our context, is valid for modes 
leaving the horizon towards the end of the anisotropic phase. 

*This is in contrast to the case where the anisotropy is present during the evolution of perturbations at any time in the 
radiation through to the present epoch. In this case the Einstein-Boltzmann system must be modified to account for the 
anisotropic evolution of the modes [241 1251 126] . 
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inverse Hubble rate due to the potential energy of (/> . For this reason, in order to have a sizeable effect, 
we need to start from a significant anisotropy. The model then admits two possible background solutions 
at early times, according to whether the expansion rate of the scale factor a is greater or smaller than the 
one of the scale factor b. In the first case (which we call positive branch), a oc t at asymptotically early 
times, while b is constant. In the second case (which we call negative branch), the anisotropic direction is 
initially contracting, a oc t~^/^ , while the other two directions expand as 6 oc t^/^ . In both cases, the metric 
is singular at t = , and the early time evolution is of the Kasner [28] type. 

In Section [3] we compute the primordial perturbations about this geometry, by extending to this case the 
Mukhanov-Sasaki [29j linearized computation and quantization procedure for the perturbations valid for a 
FRW background, [fl For clarity, we summarize at various stages the results of the isotropic computation |34| . 
and we then show how each of them extends to our case. We start by showing that also in the present case 
there are three physical modes of the perturbations. We see this with a simple counting of modes, which 
holds in the standard case, but which does not rely on the isotropy (nor homogeneity) of the background. 
For a FRW geometry, the three modes are encoded in a scalar perturbation, and in the two polarizations of a 
tensor mode (gravitational wave) . These modes are decoupled at the linearized level, due to the homogeneity 
and isotropy of the background. In our case, only one mode is decoupled, due to the residual 2d isotropy, 
while the two remaining ones are coupled to each other. As the background becomes isotropic, the latter two 
become the scalar mode, and one of the two tensor polarizations. To extract the physical modes, we choose 
a non-conventional gauge that preserves all the go^ metric perturbations. Such modes are nondynamical, 
and can be readily integrated out. Since there are initially 11 perturbations (symmetric Sg^i, and 5(1)), we 
are left with the 11 — 4(gauge) — 4 (nondynamical) = 3 physical modes. Although we perform the explicit 
computation only for the Bianchi / model with a residual 2d isotropy, we emphasize that this procedure can 
be further extended to more general geometries. In the general case, one expects that all the three physical 
modes are coupled to each other at the linearized level. 

To compute the power spectrum of these modes, we need to study their early times frequency. In the 
standard case, the physical frequency is initially dominated by the physical momentum, and it changes 
only very slowly (adiabatically) due to the expansion of the universe. In conformal time r/, related to the 
physical time thy dt = adr] the (conformal) frequency is actually constant at early times. Therefore, one 
can consistently start from an adiabatic vacuum in the asymptotic past (this, in turn, leads to a nearly scale 
invariant primordial power spectrum). The situation is more complicated for the two Bianchi / backgrounds 
that we have described above; first of all there is some ambiguity in the choice of conformal time (since 
one may use different powers of the two scale factors in the definition). This has lead us to study the 
condition under which two different conformal times can be employed, and result in the same prescription 
for the initial conditions. We find that, in the negative branch, the frequencies of two coupled modes are 
not adiabatically evolving at early times; it is possible to find a time variable in which the frequency of 
the decoupled mode is adiabatically evolving; however, the mode is tachyonic, leading to an instability and 
to a nonlinear regime that prevents us from computing firm predictions for the CMB anisotropics in it. 
For the positive branch, (o oc t, 6 — > const.), all the three modes are decoupled at asymptotically early 
times, and their three frequencies are adiabatically evolving in the conformal time variable r/, defined as 
dt = aV^ drj (where 7 any constant; for simplicity we then choose 7 = in the computation). This allows us 
to consistently start from an adiabatic vacuum. The power spectrum that we obtained approaches isotropy 
for the modes which exit the horizon right after the universe becomes isotropic. At larger scales, however, 
the spectrum exhibits an angular dependence, and it actually diverges for momenta that are aligned along 
the directions which are not expanding at asymptotically early times. As a consequence, a computation 
of the temperature anisotropics requires that we resolve this singularity. In the concluding Section [3 we 

^Perturbations in Bianchi models have been previously studied in [3D (through the formalism of [31] ') and in [3^ (through 
the formalism of !33.)- These analyses do not compute the canonical variables of the actions of the perturbations, as instead 
done in ^9] for the FRW case, which, as we shall see, is a necessary step to compute the initial conditions for these modes. 
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discuss several ways in which this can be done (at the very least, this will occur due to nonlinear effects, 
since the singularity of the modes indicates a breakdown of the linearized computation). However, all of 
them require additional input in the computation. We therefore conclude that the Bianchi / model we have 
studied is too simple to provide firm predictions for the anisotropy, at least in the assumption of initial slow 
roll of the inflaton that we have made in this analysis. 

The paper is concluded by some Appendixes, where we outline several intermediate steps of the compu- 
tations. 



2 The aim covariance 

In this Section we compute the statistical correlation between different multipoles of the CMB temperature 
anisotropics, relaxing the usual assumption of statistical isotropy. The starting point is the power spectrum 

(7^ ik)n*{k')) ^ ^^5(k - k')P(k) (2) 

of the primordial comoving curvature perturbations (we assume that the perturbations are Gaussian). If 
the background is isotropic, the fluctuations are statistically isotropic, and the power spectrum P depends 
only on the magnitude of k . In general, it will depend also on its direction. 

In the following Sections, we restrict our attention to a Bianchi / cosmology. The present computation 
is instead valid for more general anisotropic backgrounds, with the assumptions that (i) there is a residual 
isotropy of two of the three spatial dimensions, and (ii) this anisotropy is damped away during inflation. 
It is straightforward to generalize the calculation if condition (i) is relaxed. Under these assumptions, the 
primordial power spectrum depends on the magnitude of k , and on the angle between k and the anisotropic 
direction. We denote by ^ the cosine of this angle. 

We denote by (^T(p, r/o, xq) the temperature perturbation in the direction p, as measured by an observer 
at position xq and at conformal time r]Q . It is customary to decompose the temperature perturbations that 
we observe into spherical harmonics: 

'^^(p, 7?o, Xo) = ^ aijn Yim (p) (3) 

Since we are assuming that the perturbations are Gaussian, their statistical properties are still encoded in 
the second order correlations 

C«'mm' = (a^m^^'m') • (4) 

The correlations however are not diagonal as in the isotropic case. 
In the linear regime, the temperature anisotropy ([3]) is given by 

5T(p,7?o,xo) = y ^7^(k,r/i)A(A;,k•p,r/o)e''^•''^ (5) 

where the transfer function A(/c, k • p) describes the change in amplitude of the radiation perturbation from 
an initial time T]i , which can be taken deep in the radiation dominated era, to today. Since the background 
is isotropic soon after the onset of inflation, the transfer function is also isotropic (in real space). As, a 
consequence A depends only on the magnitude of the wavenumber A; = |k| of the mode, and on the angle 
between k and the line of sight at which the photon is being observed. It is convenient to expand it in a 
basis of Legendre polynomials: 

A(A:, k • p, r/o) = Y.^-i)\n + \)Pi{k ■ p)A,(fc, %). (6) 



5 



We rotate our coordinate system such that the anisotropic direction hes on the z— axis. Starting from 
the above expressions, after the algebra outhned in appendix |Xl we find 



_ 5w , M (2^ + l)(2F + l)(^-m)!(f -m)! 

c«'w - ^(-^) y (^+^)!(^/+^)! 

-A, (fc, r?o) A,, {k, r?o) j ^ d^Pr (0 (0 (A^, , (7) 

where the are associated Legendre polynomials. The above expression conveys all the statistical infor- 
mation present in the CMB in the axisymmetric anisotropic case. When constraining anisotropic models 
the full covariance should be compared to the data. 

Some general properties of the correlations follow from Eq. [71 

(i) The correlations reduce to the standard result in the isotropic case, P (/c, ^) = P {k) 



4w = l^ll'^mm' J y P (k) Ae {k, (8) 

in which different multipoles are uncorrelated. In any realistic model, the power spectrum becomes isotropic 
at small scales (large k), since such modes exit the horizon when the background is isotropic . Small scales 
correspond to large I (mathematically, this is due to the transfer functions, which are peaked at increasingly 
large k as / increases). Therefore, off-diagonal correlators become progressively smaller as / and /' increase, 
and the result ([7]) reduces to the isotropic one. The scale at which this happens is strongly sensitive to the 
amount of inflation which takes place after the universe has become isotropic. If this stage is too prolonged, 
the power spectrum is anisotropic only at too large scales to be observed today, and the isotropic result ^ 
is recovered for all multipoles. 

(ii) The residual 2d isotropy of the background restricts the number of non vanishing correlators. The 
restriction is mostly manifest in the coordinate system considered here, where the anisotropic direction 
coincides with the z— axis, and the multipoles are uncorrelated unless m = m! . For a general orientation 
of the anisotropic direction we find a greater number of non vanishing correlators. However (obviously) the 
number of linearly independent correlators remains the same. 

(iii) Bianchi / cosmologies have planar reflection symmetries (x — > — x), and, as a consequence, P{k, — ^) = 
P{k, ^). In this case, the correlation ^ also vanishes whenever the difference between / and is odd (math- 
ematically, this follows from the parity properties of the associated Legendre polynomials). 

Naturally, the orientation of the anisotropic direction is not expected to be correlated with any local 
preferred direction; for any coordinate system, the orientation can be given in terms of the Euler angles 
(a,/3, 7); these are extra free parameters in anisotropic models, which need to be fit to the data in addition 
to the usual physical parameters V determining the radiation transfer function and the spectral parameters 
S determining the primordial power spectrum. 

In principle, the anisotropic models can be fit to the observed correlations in the data by evaluating a 
likelihood of the form 

L[a|C(T',cS,a,/5,7)] = ^^— expf-iat.C-i-a), (9) 



where A'^ is the number of modes in the expansion, a is the vector of observed spherical harmonic coefficients 
and C is the model correlation (Eq. [7]) between i,m and ,m' pairs. The orientation of the anisotropic 
direction can be rotated using the spherical harmonic rotation operators I?^^/(a, /3, 7) fSS]. For the ax- 
isymmetric case one of the Euler angles is degenerate and can be marginalized out; this is not the case for 
general anisotropy. 

A full maximum likelihood search of ([9]) is computationally intensive since the dimensionality of the 
matrices involved scale as (£ -|- 1)^. In addition, cut-sky effects complicate the rotation of the data vector 
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which is computationally much faster than multiple rotations of the model correlations C. For this reason 
we postpone any detailed analysis of the likelihood surface and focus here on the framework required to 
compute the correlations in the axisymmetric anisotropic case. As we will see below, this requires a careful 
treatment of both the theory of perturbation evolution in an anisotropic background and the early-time 
initial conditions of the perturbations. 



3 The Background 

Here and in the remainder of the paper we focus on a Bianchi / model, with equal expansion rate in two of 
the three spatial dimensions: 



dt^ + a {tf dx^ + b {tf ( dy^ + dz 



(10) 



(contrary to the previous Section, we take the anisotropic direction to coincide with the x-axis). This model 
describes an anisotropic universe as long as the ratio between the two different scale factors evolves with time 
(since any constant ratio can be rescaled to one). The two physical quantities characterizing the geometry 
are the two expansion rates 

F, = - , Hb = \ (11) 
a 

where, as usual, dot denotes derivative with respect to physical time t . 

Since we want to study inflation in this background, we consider a scalar field (p, with the potential V. 
The action of the system is 



S 



M2 
lyj-p 



d^X, 



-gR 



d^Xy 



-g 



1 



(12) 



We decompose the scalar field as (p = (p + 5(j), namely in a background quantity (j) (that we take to be 
homogeneous, as the underlying geometry), plus perturbations 5<j) . In the present Section we discuss the 
background evolution, so we only include 4> here. The scalar field perturbations will be considered in the 
remainder of the paper, when we also perturb the background metric (IIOD . 
The background evolves according to the three independent equations 



H + 3H'' 
3H^ -h"^ 



V/Ml 
V = 0, 



+ 



1 



1 



+ v 



where 



H 



Ha + 2Hh 



h 



Hg — Hf) 

V3 



(13) 



(14) 



{Mp is the reduced Planck mass, while a prime on the potential denotes differentiation with respect to (f)). 
The first two equations are identical to the ones obtained in the isotropic case, in terms of the "average" 
expansion rate H . The third equation, appearing as a "modified Friedmann equation" , can be then used as 



• 2, where Eq.! 



'^"The first of p3|) is obtained by combining tlie three nontrivial Einstein equations (namely, 3Eq.Q — Eq.]^ — 2Eq.: 
refers to the ij Einstein equation) , while the third one is the 00 Einstein equation. The second one is the equation for the scalar 
field (any other linear combination of Eq.jJ, Eq.\, and Eq.j can be obtained from (|13p . as a consequence of the component 
Bianchi identity). 
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an algebraic equation for the difference h between the two expansion rates. By differentiating this equation, 
and by combining it with the other two, one can also find 



h{h + 3Hhj =0, (15) 

which is solved either hy h = (isotropic case) , or by a time evolving h which is decreasing by a rate set by 
the average expansion parameter H . 

We want to study the possible background solutions for this system under the assumption of initial slow 
roll of 4> . Let us first study the approximate case in which V is constant, and the inflaton is at rest (vacuum 
energy). The first and third of (I13p are then solved by 

H = —= ■= , /i = lb — — ■= , V constant (16) 

We note the presence of two distinct branches of solutions, characterized by either positive or negative values 
of h . In the following, we refer to them as the positive and the negative branch, respectively. 

We see that the metric has a singularity at t = . Close to the singularity, the system approaches a 
Kasner [28] "vacuum solution" , where the presence of the source (V) can be neglected. The line element of 
Kasner geometries is 

dsi^nev = -dt^ + t^'^dx^ + t^'^dy'^ + t^^dz^ (17) 

where the exponents satisfy the two properties a + /3 + 7 = + + 7^ = 1 . In the present case (/3 = 7), 
we can invert the relations ()14p . to find, close to the singularity. 

Ha ^ - , Hh ^ , positive branch 
1 2 

Ha ^ -— , Hb ^ — , negative branch (18) 

Integrating these rates reproduces the line element (|17p with a = l,/3 = 7 = 0in the positive branch, and 
a = —1/3, /3 = 7 = 2/3 in the negative branch. 

We note that the singularity is only a coordinate one for the positive branch, while it is a real one for 
the negative branch. This can be seen by computing the curvature invariants in both branches; while R is 
finite at t = in both cases, the product R^^paR^^'"^ is finite in the positive branch, while it diverges in 
the negative one. It is possible to perform a change of coordinates in the t, x plane for which the metric is 
regular in the region corresponding to t = . This would select a new time, which could be used for the 
quantization of the perturbations. However, we prefer to choose the time t for the quantization, since the 
background is constant on fixed t surfaces. As we will see, this procedure allows us to choose an adiabatic 
initial vacuum for all modes. 

We also see from (I16p that, for both branches, the system quickly reaches isotropy on a timescale Mp/\/V 
(which parametrically coincides with the inverse of the Hubble rate due to the vacuum energy), after which 



one quickly approaches the de Sitter values H = ^V/3Mp , h = 0. We note that, in the negative branch, 
the scale factor of the anisotropic direction experiences a bounce, since it contracts close to the singularity, 
while at late times it expands with the asymptotic de Sitter rate. 

Let us now return to the complete system (I13p . The only restriction that we impose is that the inflaton 
(p is initially in a slow roll regime. Since h decreases very rapidly, we need to start from a strong anisotropy 
in order to have a measurable deviation from the standard (isotropic) inflationary results. To do so, we 
choose \h\ ^ /Mp at the start. In this regime, the average rate H is much greater than in the isotropic 
case, so that, due to the higher Hubble friction, the inflaton is initially rolling even more slowly than in the 
standard case. For this reason, the results obtained above for constant V describe very accurately also the 
evolution of the complete system during the isotropization stage. 
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The degree of accuracy can be checked for any given inflaton potential, by solving the system (I13p with 
an expansion series in time. For massive chaotic inflation, with V = m?(j)'^/2 , we find 



H 



1 

37 



1+ + 



2M2 



rr? (f)Q 
4M| 



+ 



+ 



(19) 



(these expansions do not assume m <C Mp, although this is the regime of phenomenological interest). In 
figure [T] we plot the time evolution of the two expansion rates Ha and Hi,, and of the inflaton field; the 
left panel refers to the positive branch, while the right one to the negative branch; Vq denotes the initial 
potential of the inflaton. The evolutions shown (obtained by a numerical integration of the system (llSp ). are 
in perfect agreement with the above discussion. The background has an initial strong anisotropy (we choose 
/io = ±10^ m in the two branches, respectively), which is rapidly damped away, and it is then followed by 
a stage of standard isotropic inflation. Since the inflaton is nearly static in the early anisotropic stage, its 
initial value controls the duration of the isotropic inflationary expansion. The value = 16 Mp chosen in 
the figure leads to about 60 e-folds of standard infiation. 




Figure 1: Time evolution of the two expansion rates Ha (red-solid) and Hi, (green-dashed), in units of m, 
and of the infiaton (blue-dotted), in units of Mp, for a choice of chaotic infiaton potential V = m}(\?' jl. The 
left panel refers to the positive branch, while the right one to the negative branch. Vq denotes the initial 
potential of the inflaton, with c/iq = 16 Mp in this example. Notice that for the negative branch Ha starts 
from negative values, indicating that this direction is initially contracting. 



4 Cosmological Perturbations 

We now compute the perturbations about the backgrounds presented in the previous Section. 
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It is useful to recall first the standard results. In presence of a single scalar field (the inflaton) in 
a homogeneous and isotropic FRW universe one finds three physical perturbations. These three modes 
are always decoupled from each other, due to the symmetry properties of the background. The actual 
computation can be found, among other works, in the comprehensive review [33]. The two polarizations of 
the tensor mode are denoted by hx and The scalar degree of freedom is instead encoded in the so called 
Mukhanov-Sasaki variable v [29], which is in turn related to the comoving curvature perturbation TZ as 



These quantities are gauge invariant (they do not change under coordinate transformations), and TZ repre- 
sents the gravitational potential on comoving hypersurfaces where the inflaton field is homogeneous (this 
property defines these hypersurfaces). 

In the remainder of this Section we discuss the anisotropic case. We do so in three parts. In the first 
part, we show how the isotropic computation can be extended to the anisotropic background presented in 
the previous Section. We also outline the procedure to further generalize this computation to arbitrary 
anisotropic spaces. In the second part we discuss the initial conditions for the perturbations both in the 
positive and negative branch. We conclude with the computation and the discussion of the curvature power 
spectrum. 

4.1 Anisotropic background 

The number of the physical modes can be obtained by a simple counting. We start from 10 perturbations 
in the (symmetric) metric, plus 1 perturbation of the inflaton field. Of these 11 perturbations, 4 are gauge 
modes, which can be set to zero once we completely fix the freedom of general coordinate reparameterization. 
Of the remaining modes, 4 are nondynamical (namely, they enter with at most one time derivative in 
the action), and 3 physical. A convenient method to extract the 3 physical modes is to choose a set 
of conditions for the perturbations which fixes the gauge completely (thus removing the 4 gauge modes) 
without eliminating any 6gQ^ mode. These modes are nondynamical, as can be most easily seen in the ADM 
formalism [3 6 J , where the metric elements qqq and goi are written as Lagrange multipliers. We can therefore 
integrate them out of the action, as we explicitly show in Appendix IB. 31 and IB. 41 P^l We remark that this 
counting does not assume any symmetry of the 4 dimensional background. Therefore, the procedure that 
we have just outlined can be applied to more general geometries than the simple Bianchi / considered here. 
We first performed the computation in conformal time rj, defined as 



We explain the reason for this choice in the next subsection, where we also discuss how different choices of 
time can affect the quantization procedure. 

We denote the physical perturbations about this background by Hx , H^, and V, with the understanding 
that, when the background becomes isotropic, each quantity becomes the corresponding lower case pertur- 
bation introduced in the previous subsection. The explicit definition of the physical modes in terms of the 
starting metric and infiaton perturbations is presented in appendix [Bl Here we only outline the equations 
satisfied by these modes, while in the next subsection we give their initial conditions. The form of these 

^^This is analogous to what happens in electromagnetism. There, one starts from the 4 components of the vector potential. 
One degree of freedom is eliminated by fixing the U (1) gauge, while the Aq component is non dynamical (it can be fixed by 
Gauss' law). One is then left with the 2 physical polarizations of the massless photon. Clearly, one can also choose different 
gauges; for example, a gauge often chosen in the isotropic gravity computation is the longitudinal one, which removes the Sgoi 
modes. In general, however, the nondynamical modes will be linear combinations of the remaining perturbations, and it may 
be less straightforward to find them. 




(20) 
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equations is the same for both branches. The mode Hx is always decoupled from the other two, as a conse- 
quence of the residual isotropy in the y — z plane. The other two modes are coupled to each other. Denoting 
by prime a derivative with respect to conformal time, we find 

H':, +L0lHy, =0 




= 0. (22) 

The explicit expressions for the frequency elements are given in eqs. (IB.23|) and (IB.32p . 

Correspondingly, the quadratic actions for these modes are formally identical to the action of oscillators in 
Minkowski space with time evolving frequencies. This generalizes the Mukhanov- Sasaki [29] computation, 
valid in the FRW case: we can then proceed to the quantization of the modes, and provide their initial 
conditions (see the next subsection). 

As the background becomes isotropic, 6 ^ a, we find 

n" 

2 2 7 2 "" 
"^x ' "^22 ~^ ~ 

2 ,2 ^" 

u;n-^ -- , 

u;?2 ^ (23) 

where k is the comoving momentum. These are the standard evolution equations for the FRW case. All the 
modes decouple, and the equations for the two tensor polarizations become identical. 

4.2 Initial conditions 

The goal of this subsection is to set the initial conditions for the perturbations at asymptotically early initial 
times. The key quantities are the frequencies (IB.23P and (|B.32p . which we evaluate at asymptotically early 
times through the background relations (|19p . 

Again, it is useful to start by summarizing the standard isotropic computation. The physical modes 
hx, V are those that diagonalize the quadratic action for the perturbations. In conformal time, defined 
as 

ds'^ = c? (r/) \-dl^ + dx^ + dy^ + dz^l (24) 
the action for each of these modes (after Fourier transforming the spatial coordinates) is of the form 

5(2) = ^Idr^d^fcjl^f 

= k^-c?j:{ri) (25) 

where (5 denotes any of the modes, and .7-" is a time evolving function which is of the order of . The term 
q}T is exponentially small at asymptotically early times. The net result is that at early times the action of 
the mode (once written in conformal time) approaches the one of a simple harmonic oscillator in Minkowski 
spacetime, with the frequency equal to the comoving momentum k . In the Minkowski case we would have 

^Mink. = (26) 



for a choice of positive frequency modes. In the cosmological context, this is replaced by the adiabatic 
vacuum 

e 



-i drj' L 



2uj 



(27) 
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which is a solution of the evolution equation as long as the frequency is adiabatically changing, w' <C . 
The fact that k dominates w at early times, can be also restated as k/a ^ H . Going backwards in time, the 
universe is nearly exponentially contracting during inflation, and the physical momentum p = k/a becomes 
the dominant quantity for the evolution of the mode. 

Let us now discuss the anisotropic situation. We separate the comoving momentum k = ki + k2 , where 
ki denotes the component along the privileged direction x, and k2 the component in the y — z plane. We 
also denote by ki and k2 the magnitudes of the two components. The magnitude of the physical momentum 
is therefore 

/ = ^ + (28) 

For the positive branch, during the anisotropic phase, a (x t , while b is nearly constant. As we go 
backwards in time, only the x direction is "squeezing" , while the other two become frozen. Correspondingly, 
ki/a ^ k2/b at sufficiently early times, and we therefore expect that only the component of the momentum 
in that direction sets the initial condition. To be able to set an initial vacuum as in the isotropic case, we 
therefore use the scale factor a in the definition of the conformal time, cf eqs. (I2ip and ()24p . Indeed, by 
inserting the early time asymptotic behaviours (|19p into the various frequency elements of eq. (|22p . we find 

^2 _ ,.2 , ^ /_2,2^ ( ^Ix ^i2 \ _( kl + (niH^) O (mH^) 



O (mH^) kj + {mH^ 



so that ki always becomes dominant provided we go sufficiently early in timeEl Moreover, we see that 
the modes and V decouple also at asymptotically early times, so that we can set the initial conditions 
directly on them (rather than on the modes which would diagonalize the frequency matrix). Therefore, as 
long as 



/ / / 2 2 

^11 ^22 ^12 ^12 

2 ' 2 ' 2 ' 2 ' 2 

IVy, UJii LO22 ^11 ^22 



<. 1 (30) 



we can consistently start from the adiabatic vacuum ()27|) for the perturbations. This always happens, 
provided we go sufficiently early in time (the smaller ki is, the earlier we have to go in time to satisfy these 
conditions). 

In appendix IC.ll we study whether other choices of time are possible to set the initial conditions. In 
the isotropic case, the conformal time is related to the physical one by dt = adij. The most obvious 
generalization is to consider some arbitrary power of the two scale factors in the definition of the conformal 
time. We therefore discuss conformal times r related to ij by 

dri = a~2° b-^f^dr (31) 

where a and /3 are two constant parameters. We show in appendix IC.ll that this leads to an adiabatic 
evolving frequency only if a = 0. Since b is constant at asymptotically early times, we conclude that 
the time ry used here is the only possible one for the quantization, up to a trivial constant rescaling. In 
particular, this rules out the most "symmetric" choice dt = a^/^ b'^^^ dr, where the scale factors along the 
three coordinates are used with an equal weight. 

Let us now turn to the negative branch. In conformal time r], the frequency elements appearing in 
eq. (|22p are the same as for the positive branch (namely, eqs. ()B.23p and ()B.32p ). However, the early time 



'^■^The only exception is when fci is strictly zero, that is when the momentum of the mode lies in the y — z plane. This is a set 
of measure zero in the integral ((TJ for the correlators Cuimm', so we can disregard this exceptional case. Notice also that we 
used physical time to indicate the residual entries in the frequency element, while the evolution equations (|22p are in conformal 
time. 
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asymptotics are now different (due to the fact that we now take the negative solution for h in the asymptotic 
eq. dHD): 



^11 
^12 



^12 
"^22 



4 

9F 



^2 







4 fc2 



+ 0(1) 



(32) 



where we recall that in the negative branch a oc t~^/^, b oc t^^^ at asymptotically early times. 

The situation is now significantly worse than for the positive branch. Firstly, the frequencies are not 
adiabatically evolving at early times: 



\UJy 



4 



^11 



^^22 



^22 



(33) 



For this reason, we are not able to start from an adiabatic vacuum at early times, as we could do for the 
positive branch. Secondly, the 2d vector is tachyonic. 

As we show in appendix IC.21 it is not possible to find a conformal time in which the frequency of the 
2d scalars is initially adiabatically evolving. The situation is somewhat better for the 2d vector modes; by 
performing a time redefinition of the form ()3ip . with a = 2 (1 + /3), the frequency approaches a constant 
negative value in the asymptotic past. In this case, we may choose to start from an adiabatic vacuum (using 
the absolute value of lo'^). However, the fact that is negative makes the mode diverge. To see this, we 
note that, for a = 2 (1 + /?), the transformation between the physical time t and the comoving time r reads 



dt = adr] 



-2a+l T -2/3 



dr = ctdr 



(34) 



where c is a positive constant. Integrating this relation, we see that the initial singularity (t = in physical 
time) occurs at r = — oo. Since approaches a negative constant as r ^ — oo, the 2d vector experiences 
a tachyonic growth which makes it diverge. More appropriately, we conclude that the 2d vector becomes 
non-linear, and the perturbative computation breaks down. Modes with different scales and directions of 
momenta are coupled at the nonlinear level; this will likely "isotropize" the effective Hubble drag felt by the 
non-linear modes. Intuitively this would lead to a feedback resulting in the suppression of the their growth. 
Nonetheless, the breakdown of the linearized computation (together with the lack of adiabaticity in the 2d 
scalar sector) prevents us from making trustable and firm predictions for the temperature anisotropics in the 
negative branch, and for this reason we do not consider this branch further in the remainder of the analysis. 



4.3 Curvature power spectrum 

We want to compute the power spectrum for the (late time) scalar mode. According to the discussion of 
the previous subsection, we restrict our attention to the positive branch. We proceed as follows. We fix 
an initial value for the inflaton ((/>o = 16 Mp in all the cases studied in this work, so to have about 60 e- 
foldings of isotropic inflation), and we scan over several values for the two components ki^2 of the comoving 
momentum of the mode. We choose a sufficiently early time such that the conditions ()30p are satisfied, and 
we set the initial conditions for the modes and their derivatives as given by the adiabatic solution (j27p . We 
then numerically evolve the evolution equations (j22p for the coupled {V, H^} system, until a final time for 
which the background is isotropic (so that V = v, = /i+) and the modes of interest are well outside 
the horizon. In figure [5] we show a contour plot with the value of A;^|7;p at the end of inflation. This 
quantity is the curvature power spectrum, up to an overall normalization factor. The horizontal axis in the 
figure is the magnitude k = \J ki+ k2 of the comoving momentum, in units of a reference momentum A;iso. 
The latter is defined as the comoving momentum of the modes which exit the horizon when the universe 
becomes isotropic; namely, k\so = o-H at t = tiso , which we conventionally define to be the time at which 
h = 10~^m. The vertical axis is instead the cosine of the angle between the comoving momentum and the 
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Log[o(k/kiso) 



Figure 2: Power spectrum of the comoving curvature perturbation TZ in the inflationary model V = m? </>^/2. 
The inHaton field starts with i;^ = 16 Mp, and is evolved until the moment shown in the figure, when (j) = Mp. 
Modes with k > fcjso leave horizon during the later isotropic stage of inflation. The standard result is 
recovered for these modes. 

X— axis (namely, ki = k^). In figure [3] we show instead some sections of the power spectrum at fixed values 
of 

Modes with k ^ k[so leave the horizon during the isotropic phase; at early times, their frequency changes 
adiabatically, and, as a consequence, they evolve according to the adiabatic solution ()27p . As the background 
becomes isotropic, the frequency iv becomes the standard one, so that the mode adiabatically evolves to the 
adiabatic solution that one would have also found if the background had always been isotropic. For this 
reason, we expect to recover an isotropic power spectrum (no ^ dependence) in this limit. This behaviour is 
manifest in the results shown. Larger scale modes {k ^ kiso) are instead more sensitive to the background 
evolution in the anisotropic phase, and we expect nonstandard results in this regime. Indeed, at any fixed ^ , 
the power spectrum exhibits an oscillatory behaviour for k ^ kiso, and it then sharply decreases for k ^ . 

We also observe that the power spectrum increases at low ^ . The numerical results indicate that the 
the power spectrum actually diverges as 1/^ there: 

V (k) 

P{k,^)^—^ forO<e<e(A;) (35) 
? 

The region where this happens shrinks to smaller and smaller size as k grows (namely, ^ {k) is a decreasing 
function of k in the region we have probed), in agreement with the fact that the power spectrum approaches 
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PR/27r2 




Figure 3: Sections of the power spectrum (same evolution as in the previous figure) at fix values of ^ . We 
see that the power spectrum becomes isotropic (no ^ dependence) at large momenta; at lower momenta 
it presents some oscillations, whose amplitude increases as ^ decreases (this leads to the 1/^ divergency 
mentioned in the main text). Finally, the power spectrum is suppressed as A; — > . 

the isotropic result for /c — > oo . Nonetheless, the growth of the power spectrum takes place at all values of 
k we have computed, provided ^ is sufficiently small. 

Although we do not have an analytical solution of the system ()22|] . the large power at small ^ is most 
probably due to the initial normalization \via\ oc oc ^"^/^ . In the isotropic case, modes with very small 

values of k start with a large amplitude |fin| oc , but this does not lead to a large power spectrum 

P oc since the prefactor k^ is small in this regime. For the anisotropic case, the problem arises when 

k is finite, but ^ ^ . In this regime, ki <^ k2 , so that the physical momentum is almost lying in the y — z 
plane in the later isotropic regime. However, the initial condition for the mode is set at sufficiently early 
times, for which /ci/flin ^ A;2/6in . At the initial time, the mode has the large amplitude |fin| oc , while 
the prefactor k^ in the late time power spectrum A;^|z;p equals fcf , which is not suppressed in the ^ ^ 
limit. 

The divergency in the power spectrum results in a logarithmic divergency in the Cu'mm' correlators, due 
to the fact that the associated Legendre polynomials approach a finite constant value for ^ — > in the angular 
integral of eq. ([7j), while the power spectrum diverges there. We do not expect that this singularity will 
actually occur, but it rather indicates the breakdown of the linearized computation for the metric/inflaton 
perturbations. This means that the precise value of the correlators is sensitive to the nonlinear dynamics of 
the modes, which cannot be accounted for in the present linearized computation. 

In addition to this, we note that the growth (|35|) occurs at small ki . We recall that the early time 
frequency of the modes is = /ci + O (t^) , see eq. (fT9|) . and that the metric is singular at t = . The 
adiabaticity condition, uj' juP' ^ 1 is satisfied by all modes at sufficiently early times. However, the smaller 
k\ is, the closer one needs to go to the singularity for the adiabaticity condition to hold. Suppose that 
the background solution can be trusted only from some time on. Then, we can consistently start from 
the adiabatic vacuum only for modes with sufficiently high momentum k\. The value of the modes with 
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lower momentum may be simply provided as an (arbitrary) initial condition, or may be controlled by the 
additional dynamics (possibly, by additional degrees of freedom), which may be relevant at earlier times. 

To conclude, we expect that the singularity in the power spectrum is actually absent once the nonlinear 
effects and/or the complete early dynamics of the system are taken into account. In either case, however, 
this means that additional inputs are needed to provide firm predictions for the temperature anisotropics 
in this simple model. 

5 Discussion 

The main result of this work is the extension of the computation of primordial perturbations in a FRW 
geometry to more general backgrounds. We have done this in detail, only for the simplest case of a Bianchi 
/ model with residual isotropy in two spatial directions, and with one scalar inflaton field: as in the isotropic 
case, there are three physical modes. However, in contrast to this case, two of the modes are coupled to each 
other already at the linearized level, due to the fact that the background has less symmetries. At late times, 
when the background becomes isotropic, one of these two modes becomes the scalar perturbation, while 
the other one becomes one of the two perturbations of the tensor mode. The key step in the computation 
is to completely fix the gauge for the perturbations without removing the Sgofj, entries. Such modes are 
nondynamical, and can be readily integrated out of the quadratic action of the perturbations. In this way, 
one is left with only the physical degrees of freedom. It is straightforward to extend this computation to 
a general Bianchi / model, without the residual 2d spatial isotropy. In this case, all the three modes will 
be coupled to each other at the linearized level. Moreover, this procedure can be performed also for more 
general backgrounds, or for cases in which more fields are present. 

Such anisotropic backgrounds have a potentially very interesting phenomenology; for instance the early 
time coupling between the scalar and tensor modes produces an off diagonal correlation between the CMB 
modes (eq. [7]), and, possibly a nonstandard tensor-to-scalar ratio. Moreover the two polarizations of the 
tensor modes behave differently at early times, and acquire different values. Therefore, a future detection 
of the tensor perturbations (for instance, through the B polarization modes) can provide extremely useful 
information to test this possibility. 

The computation that we have outlined is a necessary first step for such studies, since it must be done 
in order to find the physical modes, and to fix their initial conditions. In standard isotropic inflation, 
the frequency of the modes is controlled by their physical momentum p at early times, which is varying 
adiabatically due to the expansion of the universe, m' /lo"^ <^ 1 . It is customary to start from the adiabatic 
vacuum (the immediate generalization of the de Sitter Bunch-Davies [37] vacuum), which is a good solution 
of the equations of motion as long as the adiabaticity condition holds, lj Due to this fact, the final result (a 
nearly scale invariant spectrum of the perturbations, in perfect agreement with observations) is insensitive 
to the "initial time" at which the initial vacuum is set, provided that the mode is deeply inside the horizon 
{p » H) at this time. 

Models which go beyond the standard FRW geometry have additional parameters. To have some predic- 
tive power, one should hope that their phenomenology depends on as few inputs as possible. For instance, if 
we find that also in this case the frequency of the physical perturbations varies adiabatically at early times, 
then the problem of setting the initial vacuum is no worse than for standard inflation. We have seen that 
this is not the case for the simple Bianchi / backgrounds that we have studied here. The main reason is that 
the anisotropy parameter h (the difference of the two expansion rates) decreases very rapidly at the onset 
of inflation. Clearly, this geometry is continuously connected to FRW, and one can take h arbitrary small 
at the initial time. In this case, however, one recovers standard FRW inflation plus negligible corrections. 

^^The choice of the initial vacuum is however object of intense debate, see for instance mostly due to the fact that the 
physical momentum is trans-Planckian close to the initial singularity, and unknown UV physics may leave some imprint in the 
mode. 
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If one hopes instead to find sizeable deviations from the isotropic case, the initial value of h must be large. 
In this work we have studied the early evolution of the model starting from \h\ ^ ^/V /Mp (where V is 
the potential energy of the inflaton), assuming that the inflaton is in a slow roll regime also during the 
anisotropic phase. This leads to two distinct possibilities: in one case, two dimensions become static at 
asymptotically early times; in the other case, the anisotropic direction is initially contracting. In the second 
case, the frequency of two of the modes is not adiabatically evolving at early times, while for the other one 
is actually tachyonic. In the first case, we do find an initial adiabatic regime; however, the amplitude of 
the modes increases and eventually diverges as their momentum is more and more aligned towards the non 
expanding directions. 

As we have mentioned, this divergence is merely an indication of the breakdown of the linear approxi- 
mation taken in our treatment. It is probable that nonlinear dynamics will stop the growth of the modes at 
this stage. However, even if this is the case, the predictive power of the model is reduced greatly, since the 
much simpler linearized treatment is invalid. There may be other physical reasons why the run-off growth 
in the mode may not occur. A better behaving vacuum may result if all the directions are expanding at 
asymptotically early times. For a Bianchi / geometry, even dropping the assumption of residual 2d isotropy, 
this cannot be achieved if the infiaton is initially in a slow roll regime (this can be understood from the 
00 Einstein equation; close to the initial singularity the three expansion rates cannot be all positive, if the 
energy density of the infiaton remains finite). Alternatively, it is possible that the Bianchi / geometry is 
only the final stage of the anisotropic phase, and that additional degrees of freedom are relevant at earlier 
times. Both these possibilities require additional input with respect to the minimal set-up studied here. 

One may also study completely different anisotropic models from the start. For instance, it may be very 
interesting to study models which evade Wald's theorem on the isotropization of Bianchi geometries [TB]. 
This can happen if one modifies Einstein gravity, for instance through the Kalb-Ramond action of string 
theory |39j . or through quadratic curvature invariants [ID]. In such case, one can find attractor solutions 
characterized by anisotropic infiationary expansion. This may avoid the run-off growth problem found in 
the model we have studied, since the difference h between the expansion rates decreases slower, or not at all, 
in these models, and therefore does not need to be extremely large in the asymptotic past. The challenge 
for such models is to find solutions for which the anisotropy decays at late times, or, alternatively, can be 
kept at a very small but controllable level. This condition itself might lead to smaller, subdominant effects 
in the data which are yet to be observed. 

We have taken a top-down approach in suggesting a model-oriented origin for the broken isotropy ob- 
served in the CMB. We have examined a particular form of anisotropy in the inflating universe but the 
treatment presented here can be readily generalized to other forms of anisotropy. If the intriguing CMB 
anomalies are verified to be of cosmological origin through more accurate, future observations, including the 
polarization data, then they may well indicate a departure from the simplest picture of infiation. In that 
case models such as the one presented here offer a minimal modification of the infiationary paradigm which 
could predict the anomalous correlations required to explain the the data. 

Note added: While this work was being completed, we became aware of the work |41j . which also 
develops a formalism for the computation of anisotropies in Bianchi I models. Although the approach is 
rather different, and a different conformal time is used in [41J, we have verified that the evolution for the 
canonical modes obtained in [5l] agree with ours in the limit in which the background has a 2d residual 
isotropy. The work [H] does not study the background solutions for the model at the extent done here, and 
it therefore does not discuss the initial conditions for the perturbations, nor the resulting power spectrum 
and CMB temperature anisotropies. We also note that some of the results for the perturbations in the 
anisotropic model discussed here had also been summarized in [i2] . 
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Appendixes 



A Explicit computation of C^'mm' 

We want to prove the result d?]) appearing in the main text. A full treatment of the line of sight calculation 
of CMB anisotropies can be found in [31]. Inverting the relation ([3|), the expression (jJl) becomes 

Cu'mm' = {aim^H'm') = j d^pdO^p {5T {p , rjo , xq) 5T {p , Tjo , xq)) Ylj^{p)Y^,^,{p), (A.l) 

The temperature anisotropy is related to the primordial curvature perturbation as in ([5]). The statistical 
correlation is therefore encoded in the primordial power spectrum ([2]) of the curvature perturbation. This 
gives 

j^-p{k) J dnpdnp,A{k,k-p,m)^^k,k-p',r^o)Y;^{p)Y,,^,{p'), (a.2) 

Note that the directional dependence of the radiation transfer functions A(A;,k • p,r/o) at conformal time 
today r/o is restricted to the cosine of the angle subtended by the plane wave unit vector k and the unit 
vector in the direction on the sky p. Since the universe is already isotropic when the initial conditions for 
the Einstein-Boltzmann system are set, the standard line of sight calculation can be used in the calculation 
of the multipole expanded transfer functions A£{k,7]Q), which depend only on the magnitude of the wave 
vector. 

We expand the transfer functions as in ([6]), and we decompose the Legendre polynomials entering in ([6]) 
into spherical harmonics 

A-jr 

^L(k •?) = -—— YIMYlm{p). (A.3) 

+ M=-L 

so that the correlator can be rewritten as 

C«w = / y^^P(k) (-0''-'''AL(A:,%)Al^,(A:,%)yLM(k)lL'Af'(k) 

X J dnpdnp,Yl^{p)Y,,^,ip')YL,,{p)Yl'M'{p')- (A.4) 

It is then straightforward to compute the angular integrals, given the orthonormality of the spherical har- 
monics. This gives 

Cii'mm' = i-iY-'' J ^^P(k)A,(A;,%)A,^(A:,r?o)r/™(k)r^w(k). (A.5) 

Note that this expression is valid for any form of anisotropic initial curvature power spectrum if the universe 
is isotropized by the time of reheating. 

By assumption, the power spectrum in the model considered here is axially symmetric around the 
anisotropic z— axis. We therefore use polar coordinates in momentum space, where 9 is the angle between k 
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and the z— axis. Due to the residual symmetry, the integral over the second angular coordinate (p is trivial. 
More specifically, we expand 

where = cos 9 , and we use the fact that 

fj d^Y,^{eA)Yl,„AOA) = 2vrW^y^^^^S^^||E^E^^ (A.7) 

to obtain eq. ([7]) of the main text. If the power spectrum is isotropic, the first angular integral is also carried 
out trivially and using the property 

we recover the standard expression ([8]). 



B Explicit computation of the perturbations 

In this appendix we perform the explicit computation of the perturbations about the background geome- 
try (I2ip . Since this is a linearized computation, we ignore the nonlinear coupling between different modes. 
We can therefore fix a comoving momentum k, and study the evolutions of the modes having that mo- 
mentum. The computation is exhaustive as long as we can solve the problem for any arbitrary value of k. 
Due to the residual symmetry in the y — z plane, we can fix fc^ = without any loss of generality. More 
appropriately, we denote by ki the component of the momentum along the anisotropic x direction, and by 
k2 the component in the orthogonal plane. We then choose the y — z coordinates such that the orthogonal 
momentum lies on the y direction. In coordinate space, this amounts to considering modes which depend 
on X and y only. The most general perturbations of the metric can be then written as 

/ -a2(l + 2$) adix 0826 

V 

In addition, there is the perturbation of the inflaton field if = (p + Sep . 

In the isotropic computation, the modes are classified into tensor, vector, or scalar modes, according to 
how they transform under 3d spatial rotations. This is particularly useful, since, due to the isotropy of the 
background, perturbations belonging to different groups are decoupled at the linearized level [HI [S]. A 
3d scalar has one degree of freedom. A 3d vector has one spatial index; however, due to the transversality 
condition {diVi = 0), it has two degrees of freedom. A 3d tensor has two spatial indices, and it is required to 
be symmetric, transverse {dihij = 0), and traceless. Therefore, it also has 2 degrees of freedom. Altogether, 
the 10 + 1 initial perturbations (of the metric and of the inflaton field), are classified into 5 scalar, 2 vector, 
and 1 tensor modes. El 

^'*In our computation we disregard the cases where the x and y dependence are triviaL Namely, we always assume that both 
fci and k2 are diflterent from zero. These two cases are of zero measure in the integral ([7| and can therefore be disregarded. 

^^We stress that this is only the first step in the computation, which is dictated by mathematical convenience. Not all of these 
modes are physical; altogether, 4 degrees of freedom are gauge modes, and 4 other ones are nondynamical. See the discussion 
in Subsection 14.11 



b^d2Es 
62(1-2S) / 



(B.9) 
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Here we do the same, but we classify the perturbations according to how they transform under 2d 
rotations (since modes wih be decoupled only if they transform differently under 2d rotations in the y — z 
plane, due to the restricted isotropy of the background). A 2d scalar has one degree of freedom. A 2d vector 
has 2 spatial indices, and is required to be transverse; therefore, it has one degree of freedom. Contrary 
to the previous case, 2d tensors instead do not exist, since there are no degrees of freedom left, once we 
impose the mode to be symmetric, transverse, and traceless. Therefore, the 10 + 1 initial perturbations are 
classified into 8 scalar and 3 vector modes. 

We conclude this discussion with a note on our notation. A 2d vector in the y — z plane is in general 
indicated with Vi , where i = 2,3. The transversality condition is then d2V2 + d^V^ = . However, as 
mentioned before eq. ()B.9p . we can impose that the modes do not depend on the third coordinate, without 
any loss of generality. Therefore, we can set V2 = , and the degree of freedom of the 2d vector is explicitly 
encoded into its component V3 . This justifies the choice of the 2d vector modes B^, B^, appearing in the 
metric ()B.9p . The parameterization of the 2d scalars entering in ()B.9p has been chosen for computational 
convenience. 



B.l Gauge choice 

We can then proceed by either constructing gauge invariant perturbations, or by choosing a gauge which 
completely removes the gauge freedom (in the standard case, this is typically done by the choice of the 
longitudinal gauge). The two procedure are completely equivalent, due to gauge invariance. We choose the 
second option, by setting 

Sgi2 = Sg23 = Sg22 = ^933 = (B.IO) 

Unlike the more conventional gauge choices, we chose not to remove any perturbation entering in the go^ 
metric elements. As we mentioned in Subsection 14.11 these modes are nondynamical [36j. Therefore, with 
this gauge choice, we know from the start that the nondynamical modes of (IB.OP are X) B, and B3 . This 
modes will be then integrated out from the action of the perturbations (see below). 

It is straightforward, although tedious, to show that this choice (i) can be made, and (ii) completely 
fixes the gauge freedom. This is true both for the isotropic (6 = a) and anisotropic background. Under the 
general coordinate transformation ^ + S,^, the metric transforms as 

^ Sg^^u - - - (B.ll) 

We parametrize ^'^ = (C'^, c?2^, ^3) , so that ^3 is the 2d vector degree of freedom, while the other 
three degrees of freedom are 2d scalars. Eq. (jB.lip then explicitly gives 
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B3 



B3 
E3 



6, 



(B.12) 



(for completeness, in the last line we have indicated how the perturbation of the inflaton transforms). Our 
gauge choice corresponds to B = T, = E = E3 = 0. It is easy to see that this is possible for one and only 
one choice of the parameters . 

B.2 Comparison with the standard modes in the isotropic Hmit (6 = a) 

It is instructive to see how the conventional modes (tensor and scalar modes of f33]) rewrite in this gauge, 
once the background becomes isotropic. This will help us interpreting our results. To do so, we start from 
the conventional mode parameterization, without 3d vectors (which are not supported during inflation) and 
with the scalar modes written in the longitudinal gauge. In the momentum space (where the comparison is 
easier) we have: 
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where k"^ = k'l + ki- 
ln terms of the general decomposition (IB. Oh (written in the isotropic limit h = a) the expression (jB.lSp 
rewrites 
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It is straightforward to write down how the final (and properly normalized) gauge invariant perturbations 
f , and /ix appear in the longitudinal gauge 
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where 6(p is the perturbation of the inflaton in the longitudinal gauge. 

We then transform from the longitudinal gauge to our gauge (jB.lOh . Namely, we start from the perturba- 



tions (|B.9p in the longitudinal gauge (|B.14p . and perform an infinitesimal change of coordinates ; 



such that the transformed modes satisfy the conditions that specify our gauge, B = T, = E = E3 = (the 
explicit transformations of the modes are given by eqs. ()B.12p in the isotropic limit b = a, b' = a'). This 
can be achieved for a unique choice of the infinitesimal transformation parameters ^'^ : 
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Starting from the longitudinal gauge (|B.14p . and performing the infinitesimal transformations ()B.12p 
with the parameters ()B.16p . the dynamical modes (5(/>, ^, B3 in our gauge become 
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We invert these relations to express the modes in the longitudinal gauge in terms of our modes. We then 
insert the resulting expressions into eqs. (|RT5]) : in this way, we find how our modes combine into the gauge 
invariant combinations v, hx- 



+ 



k'2 at 



. Mp ki k2 ~ 
-la — ■= — ; — B-x 
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B.3 2d vectors 



We now compute the quadratic action for the 2d vectors. We insert the decomposition (IB. Oh into the 
action p^ . ignoring the 2d scalars. We also set E3 = 0, according to the gauge choice (jB.lOp . and we 
expand the action up to second order in the remaining 2d modes and B-^ . This gives 
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where the boundary terms are 
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{Ha = d/a and i^b = h/h are the Hubble rates with respect to physical time, see eqs. (fTTI) ). Next, we 
Fourier transform along the spatial coordinates. We can see from the action ()B.19p that the mode -B3 is 
nondynamical, and it can be integrated out (as we discussed in 14.11 the modes 5gQ^ are nondynamical) . 
Extremizing the Fourier transformed action with respect to B^ gives 



B, = ?iB' 



p. 
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(B.21) 



where we have used for shortness the components of the physical momentum pi = ki/a and p2 = k2/b, and 
the total magnitude p 



I Pi +P2- insert this value for B^ back in the Fourier transformed action. The 
resulting action depends only on B-^ and its first derivative. We can cast it in the form 
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where 
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and where the newly introduced boundary term is 



drj d?k ■ 



-Hb + {Ha — Hb) 
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Namely, the action ()B.22p . disregarding the boundary terms, has been written in terms of the canonically 
normalized combination Hx- This mode coincides with the standard mode /ix in the isotropic limit, as can 
be seen by comparing its definition with the third of (IB.ISP when b = a (the modes actually differ by an 
unphysical phase). Its equation of motion, appearing in eq. ([22]) of the main text, is immediately obtained 
from (IE22|). 



B.4 2d scalars 



We now compute the quadratic action for the 2d scalars. The computation proceeds analogously to the 
one of previous subsection, although it is more involved due to the presence of a larger number of modes. 
We start from the metric (jB.9p . in the gauge (|B.10p . and we now ignore the 2d vectors. We also include 
the perturbation of the inflaton field in the present computation. We insert this in the action ()12p and we 
expand it to quadratic order in the perturbations. We find 
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where we stress that (j) is the derivative of the background inflaton with respect to physical time, and where 
the boundary terms are 
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We now find that the modes <I>, x, and 5 are nondynamical (they are the Sqq^ modes), while the two 
remaining modes, ^ and 6(1) are dynamical. We can then integrate out the nondynamical perturbations 
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analogously to what we did for the 2d vectors, and obtain an action in terms of the two dynamical modes 
only. 

After Fourier transforming the spatial component, extremizing the action wrt gives 



Inserting this back into the Fourier transformed action gives 



(B.27) 
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We extremize this action wrt B*. The resulting equation of motion can be written as an equation for $ 
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When we insert this back into (|B.28p . also B disappears from the action. The resulting action can be finally 
cast in the form 
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where the canonically normalized combinations are 
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We see by comparing with the first two of (jB.lSp that these modes reduce to the analogous (lower case) 
ones in the isotropic limit (up to an unphysical phase for Hj^ ). 

The term entering in the action is the square frequency matrix 
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The equation of motion for the 2d scalars, appearing in eq. ()22p of the main text, is immediately obtained 
from (fR3(l . 

The second boundary term is 
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Dependence of the frequency, and of the adiabaticity condition on 
the choice of conformal time 



In Appendix [Bl we compute the frequencies of the perturbations in conformal time r] , defined in eq. 
Here, we study how the action of the modes changes when we use another conformal time r . 

Rather than studying the 2d vector and scalar sector separately, we can give a more compact presentation 
by considering N canonically normalized fields Xi (after Fourier transforming the spatial coordinates), with 
the action 



2 



dr] d^k 



dX^ dX 
dr] dr] 



X- 





( ^1 \ 


, x^ 


X2 







{CM) 



where Q,^ is a hermitian matrix (in the present application, 
sector, respectively). 



= 1 and 2 for the 2d vector and the 2d scalar 
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Under the time redefinition 
this action becomes 



S 



drd^k 



drj = /-2 (r) dr 



(C.35) 



(C.36) 



With this new time variable, the fields Xi are no longer canonical. We can however introduce new canonical 
fields 

Yi = fXi (C.37) 



In terms of these fields, the action reads 
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dr d-^k 
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where the new square frequency matrix is 



while the boundary term is 



B 



1 



III 

f dr^ 



- I drd'k^f^^Y^Y 



(C.39) 



(C.40) 



dr \f dr 

In principle, the two actions (jC.34p and (|C.38p lead to two different quantization procedures for the 
fields, and to two different initial conditions for the cosmological application. We want to determine under 
which conditions these two procedures lead to the same physical result. Ef] Let us now assume that the 
off-diagonal terms of both 0^ and can be neglected at asymptotically early times (as this is the case for 
the modes studied in the present work, see subsection 14. 2p . As discussed after eq. ([25]) . if the frequency 0^ 
is adiabatically evolving, we can set the initial conditions (up to an irrelevant phase) 
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y/2n^ ' V dr] 
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where fJ^ here refers to the diagonal entry corresponding to the mode X^. 
On the other hand, if 0,r is also adiabatically evolving, we can also impose 
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The initial conditions (IC.42P can be used to set the initial conditions for Xi and dXi/dr] through the 
relations (|C.35p and (|C.37p . Proceeding in this way, we find 



(dYi_ 
V dr 



T^ dp' 



( dX, 
V dv 



f d^f 



1/4 



^1/4 



J dr 



f d'f 



-1/4 



(C.43) 



The physical perturbations will be linear combinations of Xi or, equivalently, Yi. We can equivalently compute the evolution 
of these physical modes by writing and solving the Einstein equations in terms of the fields Xi or Yi; the only difference that 
may arise from using one or the other set of fields is in the initial conditions that the two quantizations (either in terms of the 
time variable r; or r) impose. 
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If these conditions agree with those obtained by using the fields Xi and the time rj directly, namely 
eqs. ()C.4ip . then the two quantizations are equivalent. This happens only if 
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at asymptotically early times. One can verify that these conditions are equivalent to 
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(C.45) 



Let us now discuss under which circumstances these conditions are met. As we mentioned, we can 
consistently set both (IC.4ip and (IC.42P only if both 0^ and 0,r are adiabatically evolving (each with respect 
to its own time). Eq. (jC.39p can be differentiated to give 
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Barring accidental cancellations between the different terms, both frequencies are adiabatically evolving 
only if 

1 d-f 



f dr]^ 
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These conditions also imply that the two quantization procedures are equivalent, cf. eqs. (IC.45p . 

To summarize, we can consistently use the two different time variables r/ and r, provided the relation 
among them varies sufficiently slowly. If conditions ()C.47p are met, and if the frequency is adiabatically 
evolving, then also the frequency $7,- is adiabatically evolving in its own time. These two procedures lead 
to equivalent initial conditions for the physical modes. On the other hand, we see that not all choices for 
conformal time are appropriate for the quantization. 

We proceed by discussing what these findings imply for the Bianchi / backgrounds studied in the main 
text. For sake of clarity, we study the positive and negative branches separately in the two following sections. 



C.l Positive branch 



We show in Section [4.21 that, for the positive branch, one can quantize the perturbations in the conformal 
time r/, defined in eq. (j2ip . Here, we want to investigate whether some other time variable r can be used. 
For definiteness, let us restrict our attention to time redefinitions of the form 



dii = (r) dT 
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where a and h are the two scale factors, and a and (3 two constant parameters (this is the most immediate 
generalization of what is done in the isotropic case). Eqs. (|C.47p give the three conditions that / needs to 
satisfy for the two choices of time to provide the same initial conditions. We perform the computation using 
the time r] , and we then present the early time expansion in physical time: 
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(t) Hi it) [l + O (m^ t^) 
{t)Hl{t) [l + Opt^) 



(C.49) 

Namely, these three conditions coincide. Ell From eq. ()29p we have = kf+O {rn? t"^) for all the three modes, 
where ki is the longitudinal (comoving) momentum of the mode in the anisotropic x— direction. Therefore, 
we must have \a\a (t) Ha (t) <^ ki at asymptotically early times. The product a Ha is nearly constant all 
throughout the anisotropic sage, and very close to the value kiso defined in Section Therefore, the three 
conditions (jC.49p become 

ki > |q| kiso (C.50) 

The modes with ki ^ fcjso are those leaving the horizon after the universe has become isotropic, and 
we recover the standard result for them (see Section [4.3p . The modes with smaller longitudinal momentum 
are those sensitive to the evolution of the universe during the anisotropic stage, and therefore are those of 
real interest for the present analysis. We can only quantize them using a time variable for which a = 
in eq. (IC.48p . We remark that b is constant at asymptotically early times. Therefore, if we use the scale 
factors in the definition of the conformal time, the time r] used in the main text is the only possible one, up 
to a trivial constant rescaling. 

C.2 Negative branch 

Let us start from the early time frequencies for the modes canonically normalized with respect to conformal 
time r]. As shown in eq. ()32p we find 
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where, in the 2d scalar case, fi^ refers to the dominant diagonal entries of the frequency matrix (which 
are identical). We recall that a oc t~^^^, b oc. i^l"^ at asymptotically early times (as always, although the 
frequency refers to the conformal time r/, we show the expansion in terms of physical time t). As discussed 
in the main text, these frequencies are not adiabatically evolving at early times; moreover the 2d vector 
mode is tachyonic. Let us discuss whether the situation improves if we use a different time variable r; as 
for the positive branch, we consider time redefinitions of the form (IC.48p . 

Let us discuss the 2d vector mode first. Evaluating the two equations ()C.39p and (jC.46p we find 
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V(a- 2^ + 1) (a -2/3 -5) 

We can obtain an adiabatic initial evolution of the frequency if we choose the parameters of the transfor- 
mation such that Q = 2 (1 + /3). However, in this case, the frequency ^Ir becomes 

Q-8/3-65-4/3 
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(C.53) 



^''To be rigorous, we showed that the three conditions (|C.47|I are sufficient to ensure that the two quantizations are equivalent, 
but not necessary, since there may be cancellations in eq. (|C.46|I : however, the specific form of (|C.49|I shows that these conditions 
are also necessary, once we restrict our attention to time redefinitions of the type (IC.48[) . 
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By using the early time asymptotics for the scale factors, we can see that the frequency ()C.53|) approaches 
a constant negative value close to the initial singularity. This leads to a tachyonic growth of the mode, and 
to the breakdown of the linearized computation (see the main text). 

Let us finally discuss the 2d scalar sector. For generic values of a and /3 we now find 
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For a 7^ 2 (1 + /?), the second expression evaluates to 2 at asymptotically early time. For a = 2 (1 + /?), 
the dominant term vanishes both at the numerator and at the denominator, and the adiabaticity condition 
approaches 1/ (3tp2)- This quantity diverges as as the time approaches the initial singularity at t = . 
Therefore, no time variable r chosen as in (IC.48P leads to an initial adiabatic evolution for the 2d scalars. 
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